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Abstract 

Belle Collaboration has recently observed a new state, the X(4160), in the process of 
double charm production e+e~ J / il^ + X (4160) followed by X(4160) D*D*. We dis- 
cuss possible interpretations for the X(4160) based on the NRQCD calculations and the 
potential model estimates for the charmonium spectrum. We first focus on the D-wave 
spin-singlet 2 ^ charmonium ^D2{2D), which is estimated to have a small production 
rate of about 5% of that for e~^e~ — > J/ip + ric{lS), and therefore is incompatible with the 
observed data for X(4160). We then discuss the possibility that the X(4160) is the known 
J^'~" = 1 charmonium state ■0(4160), which can be produced via two photon fragmen- 
tation, but the production rate is much smaller than observed for e~^e~ J / ^p + X {4:160) . 
In contrast to above two possibilities, the r/c(4S') assignment is a likely one, which is sup- 
ported by the observed relatively large production rate and non-observation of DD decay 
of X(4160), but we have to understand why rye (45) has such a low mass, which deserves 
further studies. The P-wave excited state Xco(3i^) is also an interesting candidate, if 
the observed broad peak around 3.8-3.9 GeV in the recoil mass of DD against J ftp in 

+ e~ J/ip + DD is due to the Xco(2-P) state. Measurements of production angu- 
lar distributions will be helpful to distinguish between r]c{4S) and Xco(3-P) assignments. 
Production mechanisms in nonrelativistic QCD are emphasized. 



PACS numbers:13.66.Bc, 12.38.Bx, 14.40. Gx 



Using a data sample of 693 fb~^ collected around the T{4S) with the Belle detector at the 
KEKB e"'"e~ storage rings, very recently the Belle Collaboration has reported some new results 
for double charmonium production in e~^e~ annihilation at y/s = 10.6 GeVp!]. In the measured 
processes e~^e~ — > J/-ipD^*'^D^*\ a new resonance state, called the X(4160), is observed with 
a significance of 5.1 a in e^e^ J/iIjX{4160) followed by X(4160) — >■ D*D*. As a hadronic 
resonance, the X(4160) has the following mass and width[T] 

M = 415Qtfo ± 15 MeV, T = 139^11^ ± 21 MeV, (1) 

and production cross section 

a(e+e- J/0X(416O))i?^.^-. = (24.71^^38 ± 5.0) fb, (2) 

which is large and comparable to the observed cross sections for + ^ J/ip + rjd^S) and 
e"*" -|- e~ — >■ J/ip + Xco(l-P)- Although at present the data for the X(4160) are still preliminary 
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and more data are apparently needed to identify the nature of this new state, it is worthwhile 
to discuss its possible assignments, especially in view of the great potential of finding new 
particles, e.g. the ric{2S) and the X(3940)(for recent reviews on new hadrons with heavy 
quarks, see, e.g. [21 [3]) in the e"'"e~ annihilation processes at B factories. It is also interesting to 
study the production mechanisms of those charmonium or charmonium-like states in the double 
charmonium production processes in e^e~ annihilation, since the theoretical understanding for 
the double charm production is very intriguing but not totally conclusive. 

In the following, we discuss some possible interpretations for the X(4160), in connection 
with the double charmonium production problem. 

In general, in the process e^e'^ J/iIjD*D* the D*D* system can have charge parity either 
C=+ (if e~^e~ annihilated into one photon) or C=- (if e^e~ annihilated into two photons). In 
the case of C=+, the X(4160) can have J^^ = 0++, 0^+, 1"+, 2~+, 2++, while in the 
case of C=-, it will have J^'" = 1 ,2 ,1+^,2+^,... Because the two photon processes are 
relatively suppressed by an additional electromagnetic coupling constant a = 1/137, at the B 
factory energy y/s = 10.6 GeV the one photon processes usually have larger rates, and should 
therefore be considered firstly. 

I. The assignment that the X(4160) is the D-wave spin-singlet charmonium state ^D2{2D) 
is disfavored by the too small theoretical rate of production in e~^e~ annihilation. 

This 2D state has quantum numbers J^'" =2 ^. There are some arguments which could 
be in favor of this interpretation. 

First, the observed X(4160) has the same mass as the ?/'(4160) (see the Particle Physics 
Booklet [4j), which is known to be the good candidate of the D-wave spin-triplet charmonium 
state ^Di{2D) with J^^ = 1— , and 

M = 4153 ± 3 MeV, P = 103 ± 8 MeV, Pee = 0.83 ± 0.07 KeV. (3) 

The hyperfine splitting between the center of mass of ^Dj{2D) states and the ^D2{2D) is ex- 
pected to be vanishing if the short-range spin-dependent forces are due to one-gluon exchange, 
and the fine-splittings between ^Dj{2D) states should be a few tens MeV. Therefore the ob- 
served mass of the X(4160), M = 4156^20 ^ 15 MeV, could be compatible with the D-wave 
spin-singlet charmonium state ^D2{2D), and roughly speaking, is in agreement with the poten- 
tial model predictions (see, e.g. [SIEIIT]). (Note, however, that the S-D mixing effects should 
be considered due to a sizable leptonic width observed for the ■?/'(4160).) 

Second, for the J^^ = 2~~^ 2D charmonium state, the decay to DD is forbidden, while decays 
to D*D* and D*D + c.c. are allowed (also including the D*Ds mesons). Therefore, the decay 
X(4160) D*D* could be substantial. In ref.[7] the calculated decay rate of 2-+(2D) D*D* 
is nearly equal to that of 2~~^{2D) D*D + c.c. (However, the sensitivity to the model and 
parameters need to be further investigated.) 

However, the main problem for this assignment is that the D-wave state is expected to have 
a much smaller production rate than the S-wave state such as the rjc in double charmonium 
production. The double charmonium production in e~^e~ annihilation at B factories has been 
studied in the framework of nonrelativistic QCD (NRQCD) [HI |9l |10l [11] , in which the charmo- 
nium states are treated as nonrelativistic bound systems, and the production rates can be fac- 
torized into the short distance part, which can be calculated in perturbative QCD, and the long 
distance part, which can be related to the wavefunction of the charmonium. Experimentally, 
double charmonium production processes + e~ — > J/ip + ?7c(15')(?7c(2S'), Xco(l-P), -^(3940)) 
have been observed by Belle and BaBar [121 [131 [H], but the cross sections are larger than 
the leading order (LO) NRQCD calculations by almost an order of magnitude [HI [SI [H] (Note 
that the numerical results can be somewhat different when taking different parameters e.g. 
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in |H] and in [5] but the physical conclusion is the same). The next to leading order (NLO) 
QCD radiative corrections are found to be very significant to increase the cross section of 

+ e~ ^ J/ip + ?7c(15')p^. Moreover, the relativistic corrections further increase this cross 
section[16l [17]. As a result, the calculated cross section of e"*" + e~ J/i^ + Vci^S) with 
both NLO radiative and relativistic corrections in NRQCD may reach the lower bound of the 
experimental values, and could resolve the problem. 

Other approaches including the light-cone methods are also discussed in the literature to 
resolve the discrepancy between experimental data and theory for the double charmonium 
production [18]. 

It is interesting to point out that although the LO results in NRQCD for the double charmo- 
nium production cross sections in e^+e~ annihilation are much smaller than data, the predicted 
relative rates seem to be consistent with the measured values. For instance, the predicted cross 
sections for J /tp + rjdlS) and J /ip + Xco(l-P) are comparable and much larger than that for 
J/iIj + Xci(l-P) and J/iIj + Xc2(l-P)- These ratios of LO cross sections are indeed compatible 
with data. In the following, we will use the calculated LO result for e'^ + ^ J / 'ip +^ D2 and 
compare it with that for + e~ J /ip + rjc, to make some predictions. 

To the leading order in NRQCD with QED contribution included, the cross section for 
e+ + e~ ^ 7* ^ ipijiS) +^ D2{mD) process can be expressed as[TO] 
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and 



,2 40967r2(s - 16m2)3(32am2 + 96a,m2 + 3sa)^{x^ + 1) 
Mj/^^r>.\ = , (5) 



where x = cos 9, and 9 is the angle between the beam axis (pi) and the J/ip momentum (ps). 

As in[9l|10], we take following parameters: ^/s = 10.6GeV, nic = l.SGeV, ms = m4 = 2mc 
(in the nonrelativistic limit), as = 0.26, and the wave functions at the origin are taken from a 
potential model calculation (see e.g. the QCD (BT) model in Ref.[l9]): |i?is(0)P = 0.810GeV^ 
and |i?ifl(0)|2 = O.OlSGeV^ |i?2i3(0)P = 0.024GeV^. 

For the J/ip +^ -^2(1-0) (seeflD]) and J/ip +^ D2{2D) production, we then get the angular 
distributions (differential cross sections) and cross sections, which are shown in Table I, where 
9 is the angle between the incident beam and the J/ip, and the numbers with (without) square 
brackets mean the cross sections without QED (with QED) contributions. These cross sections 
are much smaller than that predicted for the J/ip + ric{lS) production [9], [10], which is also 
listed in Table I. We see that the cross section for J/ip +^ D2{2D) production is predicted to 
be only about 5% of that of J /ip + rjc, in contrast to the observed production cross section for 
the J/ip + X(4160) shown in eq.(2), which is comparable to that of J/ip + r]c [IS [T3l [T^ . 

In principle, we could also detect the ^D2{1D) charmonium, which should lie around 
3.8 GeV, in the e"^ + J/ip +^ D2{1D) process. However, the main decay modes of 

^D2{ID) should be decays to light hadrons via intermediate gluons, since the ^D2{1D) is ex- 
pected to lie below the D*D threshold. Without a dominant exclusive decay channel like D*D* 
or D*D + c.c, it will be even more difficult to detect this charmonium state especially when 
the production cross section is small. 

To sum up, although the ^D2{2D) charmonium could be a possible assignment for the 
X(4160), the predicted small production rate for + e~ ^ J /ip +^ D2{2D) makes this assign- 
ment very unlikely. Despite of the existing uncertainties in the theoretical calculation (e.g., the 
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chosen parameters, and high order corrections), this conclusion should hold, since the small 
number of 5% for the ratio of J /ip +^ D2{2D) production cross section to that of J/ip + rjc can 
not be enhanced to close to the observed value (about 1) by changing the parameters or includ- 
ing the NLO QCD corrections (Note that the NLO QCD correction to the J/ip + r]c increases 
this production rate by a factor of about 2[ 



II. The possibility that the X(4160) is the known J^'^ = 1 charmonium state '?/'(4160) 
should be ruled out. 

The '?/'(4160) is in the same mass region as the newly observed X(4160), and their widths 
are also comparable. Moreover, the ip{4160) can also decay to D*D*. However, the process 
e~^e~ J/ip + ip{4:160) can only proceed through e~^e~ annihilation into two photons due to 
the conserved charge parities. 

In fact, the two-photon process was first studied for e~^e~ — > 27* J / ip + J / ip m. ref . [20] . 
and it was found that the production rate is comparable to or even larger than that the one- 
photon process e^e' ^7* J /ip + rjc in the leading order calculation [20j. Moreover, for 
the inclusive double charm production process e~^e~ —^J/ip + cc, the two-photon process 
e~^e~ 27* — i> J/ip + cc will prevail over the one-photon process e'^e~ ^ '-/* ^ J/ip + cc 
when ^/s becomes larger than 20 GeV[21]. This is because, in these two-photon fragmentation 
processes the virtualities of the photons are only about 4m^, which is much smaller than the 
virtuality s in the one-photon process. 

However, because the ip{4160) is expected to be a D-wave {^Di{2D)) dominated charmonium 
state (with possibly some ^Si{3S) admixture), its coupling to the photon is suppressed by the 
factor |sm0_R35(O) — cog0 ^^^g -R2r'(0)P; compared with \Ris{0)\'^ for the J/tjj. Here, we have 
assumed that the '?/'(4160) is a mixture of the ^Di{2D) and ^S'i(3S') states with being the 
mixing angle: 

|V^(4040)) = |3^Si) cos (p + |2^Di) sin 0, (6) 
|V'(4160)) = -iS^Si) sin0 + \2^Di) cos0. (7) 

The above expression is only a very rough approximation, since admixtures with the charmed 
meson pairs due to coupled channel effects and with other S-wave states are all ignored. With 
this simple assumption we get leptonic decay widths for the '0(4040) and -0(4160): 



r(0;(4O4O) ^ C+C-) = ^o?e^ 



|cos0i?35(O) + s^^0i7t|^2D(O)P 



(2m 
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, JS^'^0^35(O) - COS0^^i?2^(O)|2 

r(0(416O) e+e-) = ^c?el (2m,Y • 

Using the experimental values Y^^{^[mm)) = 0.86±0.07 KeV and ree(0(416O)) = 0.83±0.07 KeV, 
and |i?35(0)|2 = 0.455 GeV^, |^2d(0)P = 0.024 GeV^, we get the mixing angle from the ratio 
of these two leptonic widths: 

= -35°, = +55°. (10) 

The mixing angle is unexpectedly large, and this is due to the observed largeness of the leptonic 
decay width of -0(4160) (almost equal to that of the -0(4040)). In fact, if we neglect the 
contribution from the 2D component of the -0(4160), we would get an estimate for the mixing 
angle that is independent of potential model parameters: ~ ±45°, which would be the 
maximum mixing. The large 3S-2D mixing is a puzzling problem in understanding the nature 
of -0(4160). Other studies like the strong decays to D*D* may be useful to clarify the 3S-2D 
mixing problem for the 0(4160) (see, e.g. discussions in [22]). 
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Despite of the above uncertainty concerning the 3S-2D mixing, we may have a quite rea- 
sonable estimate of the production cross section of e~^e~ ^ J/ip + ip (4:160), as compared with 
that of e~^e~ — >■ J/ip + J/4'- the nonrelativistic hmit, the charmonium masses are all ap- 
proximately set to be M = 2mc (i.e. all binding energies are neglected), and then we will have 
a simple relation 

a(e+e- ^ J/^ + 7/>(4160)) _ T,,{ij{mO) 

a(e+e-^J/ij + J/ij) TM/i^) ^ ' ' ^ ^ 

where the observed values T^^{J/iIj) = 5.55 ±0.14 ±0.02 and ree(V'(4160)) = 0.83±0.07 KeV 0] 
are used. This relation is obtained by the observation that in the double vector-charmonium 
production via two virtual photons in e~^e~ annihilation at ^/s = 10.6 GeV the photon frag- 
mentation is dominant (see e.g. [20l El]), in which the virtual photon converts directly into 
the vector-charmonium, the same way as the leptonic decay of the vector-charmonium. As the 
most favorable mechanism with the minimal photon-virtuality, all vector charmonium states 
(e.g. J /ip,ip {2 S),ip {4040), ip {4160), ...) are expected to be produced from the two photon frag- 
mentation in e"'"e~ annihilation at y/s = 10.6 GeV or even higher energies. 
In ref.[T3], the following upper bound is given 

a{e+e- J /i) + J /i)) x B{J/i) ^> 2 charged) < 9.1 fb, (12) 

which will imply 

a(e+e" ^ J/ip + V^(4160) x 5(^^(4160) ^> 2 charged) < 1.4 fb, (13) 

assuming B{iIj{4160) ^> 2 charged) is comparable to B{J/ip — >> 2 charged). 

This predicted cross section is much smaller than the experimental value given in eq.(2). 
Therefore, the 'il){4160) assignment for the X(4160) should be ruled out. 

III. The X(4160) could be an excited O""*" charmonium state: the ric{4S) (less likely to be 
the r/e(35)). 

As a possible candidate of the ^ state, the X(4160) can be the ric{4S) charmonium, which 
is expected to decay into D*D* and D*D + c.c, but not DD. 

Note that Belle already found a new state, the X(3940), in the process e"*" + J/i' + 

X(3940)[25j, which has a dominant decay mode into D*D (with the fraction of X(3940) decays 
with more than two charged tracks in the final state into D*D being (96^*^32 ±22)%), and a quite 
narrow width F = 39 ± 26 MeV. This result has been further confirmed by Belle (see [T|. The 
X(3940) is considered as a good candidate for the ric{3S) (for discussions see, e.g. [261 121 E]). 
The problem is the low mass of X(3940) as the ric{3S), compared with the ip{3S) candidate 
-0(4040). But this could be explained by the coupled channel effects that the coupling of ric{3S) 
to the 0"^ and 0~ charmed meson pair (in S-wave) will lower the mass of ?7c(3S')[26j. 

If we accept X(3940) as the r]c{3S), then X(4160) should be the r]c{4S) if it is a 0~+ char- 
monium. In this case, the mass difference between r]c{4S) and ric{3S) would be only 220 MeV 
. This mass difference is smaller than that predicted by the potential models with linear plus 
Coulomb potentials (see, e.g. [H [6l [H [26]). Note that the corresponding mass difference be- 
tween the ip{4S) and 4'{3S) is about 375 MeV if the ip{4S) is identified with the ■?/)(4415) and 
the ip{3S) with the ^/'(4040) as conventionally classified in the charmonium spectrum. An even 
more puzzling problem is the mass splitting between the t]c{4S) (if identified with X(4160)) 
and the ip{4S) (if identified with il^{4415)), which is as large as 255 MeV, compared with the 
mass differences between ric{lS) and J/4'{1S), rjc{2S) and ip{2S), rjc{3S) and ip{3S), which are 
only 117, 48, and 100 Mev respectively (assuming the X(3940) is identified with rjc{3S)). In 
simple potential models the mass splittings between 0~^{nS) and 1 {nS) (n=l,2,3,4,...) are 
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expected to be decreased as n increases. Although the mass spectrum can be modified by the 
coupled channel effects and S-D mixing, such a big mass difference, 255 MeV, between ?7c(45') 
and iIj{^S) is still difficult to understand, unless the assignments for excited 1 states are 
changed in some way. For instance, if the '?/'(4415) is not identified with the ip{^S) but with 
the ip{5S), as discussed in the potential model with color screening effects (see, e.g. [221 1^)- 
then the corresponding 77c(45') mass could be lowered. In this case, all higher excited states 
will be lowered in the mass spectra. But this is only a plausible resolution for the problem in 
the r]c{4:S) assignment of X(4160), other approaches apparently need to be studied. 

Could the X(4160) be the ric{3S)7 If so, what assignment will be for the X(3940). Moreover, 
if so, as the ric{3S) the mass of X(4160) would be higher than that of ipi^S), which is identified 
with '?/'(4040), by 120 MeV. The positive and large mass splitting between 0^^(35") and 1 {3S) 
seems not acceptable in charmonium spectrum. So, X(4160) can not be the r]c{3S). 

The ric{4:S) interpretation for the X(4160) is a likely one in view of the large production 
rates of ric{lS), ric{2S), and ric{3S) (if identified with the X(3940)) associated with J/ip in 
e~^e~ annihilation. For the t]c{4:S) production, the angular distribution and cross section is 
shown in Table I. Note that the normalization can be substantially enhanced with the NLO 
correction (see [TB]) but the angular distribution remains unchanged in this case. The form 
of (1 + cos^O) for the angular distribution of this assignment differs markedly from another 
interesting assignment, the 0+"^ charmonium (i.e. the Xco(3-P) state), which will be discussed 
below. 

IV. The X(4160) might be an excited 0"'""'" charmonium state: the Xco(3-P) (unlikely to be 
the Xco(2P)). 

As a possible candidate of the O"*"^ state, the X(4160) might be the Xco{3P) charmonium, 
which is expected to decay into D*D* and DD, but not D*D + c.c. The DD decay mode of 
the X(4160) has not yet been seen so far. The mass of X(4160) immediately indicates that it 
is unlikely to be the Xco(2-P) state, which is predicted to he around 3.9-4.0 GeV. The fact that 
the observed Z(3930) can be identified with the Xc2(2-P) state[4j) also implies that the Xco(2-P) 
should lie well below 4160 MeV. So, X(4160) can only be the Xco{3P) if it is a O"*"^ charmonium 
state. However, if in + e~ annihilation both J / ip + Xcoi^P) and J / ip + Xco{3P) are observed, 
why J/ip + Xco(2-P) is in the absence? In fact, according to the NRQCD calculation, the 
production cross section of J/ip + Xco(2-P) should be comparable to or even larger than that of 
J/ip + Xco(l-P) (see |TD]), because the first derivative of the wavefunction at the origin for the 
2P-state is usually larger than that for the IP-state: |-R2p(0)P > l-Rip(0)P (see, e.g.[T9]). To 
LO in NRQCD the cross sections for e+ e" ^ J/ip + Xco(l^) and e+ e" ^ J/ip + Xco(2P) 
are predicted to be 6.9 fb and 9.4 fb respectively (QED contributions are included) |10]. So, 
the experimental absence of + e~ ^ J/ip + Xco(2-P) would be hard to understand. 

However, at this point, it is interesting to notice that Belle has observed a broad peak (but 
with only 3.8 cr) around 3.8-3.9 GeV in the recoil mass of DD against J/ijj in the process 

+ e~ J/ip + D + D[l\ (it may also be seen in the 77 — > DD process). Is this the missing 
Xco(2P) state? If the bump in the 3.8-3.9 GeV region is really due to the Xcoi^P) DD decay 
, the Xco{3P) assignment for X(4160) would be favored (but the Xcoi'^P) state should be further 
examined experimentally). 

As discussed so far, two likely assignments for the X(4160) are the ?7c(45') and Xco{3P) 
charmonia. How to distinguish between them? One effectible way is to measure the angular 
distribution of the cross sections. The differential cross section in the case of Xco{3P) is shown 
in Table I (see also [lO]). Compared with (1 + cos'^6) in the case of ric{4:S), the form of (1 + 
0.252005^6*) for the Xco{3P) has a much weaker 9 dependence, and therefore the measurements 
on angular distributions can be used to test the two possible assignments. 
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V. The X(4160) is unlikely to be the excited 2+"'" or l"*""*" charmonium state, Xc2(2-P, 3P) or 
Xci(2P,3P). 

Experimentally, both e+ + e^ J / ip + Xc2{^P) and e+ + e^ J/ip + Xci{^P) have not been 
seen. This is consistent with the smallness of calculated cross sections for them. In fact to LO 
in NRQCD the cross sections for J/ip + Xc2(l-P) and J/i' + Xcii^P) are predicted to be 1.8 fb 
and 1.0 fb respectively (QED contributions are included) [IDj, which are much smaller than that 
for J / xjj + Xcoi^P) and J/ip + ric{lS). In contrast, the observed cross section for J / ip + X (4160) 
is comparable to that of J/ip + Xcoi^P) and J/ip + rjc{lS). In view of both the experimental 
non-observation and the calculated smallness of the cross sections for + e~ — > J/ip + Xc2{'^P) 
and e"*" + J /ip + Xci(lP) we conclude that the X(4160) is unlikely to be the excited 2^^^ 

or 1++ charmonium state, Xc2(2P, 3P) or Xci(2P, 3P). 

Note that the Xc2(2P, 3P) and Xci(2P, 3P) interpretations for X(4160) are disfavored by the 
experimental absence of Xci,2(lP) not only in the exclusive double charmonium production of 
Xci,2(lP) associated with J/V', but also in the inclusive prompt production of Xci,2(lP) in e+e~ 
annihilation |27j. In fact, recently Babar finds no evidence for prompt Xci,2(l-P) production 
after subtracting the contributions from prompt ^(25") production feed-down to y^i 9(1P)P7]. 
Therefore, the assignments of X(4160) as Xc2(2P, 3P) or Xci(2P, 3P) states are very unlikely. 

VI. Non-charmonium interpretations for the X(4160): glueballs, hybrids, and charmonium- 
molecules. 

As suggested in a 0"*"+ glueball associated with the J /ip may be produced with a sizable 
rate in e^e~ annihilation. However, as a O^"*" glueball, the X(4160) would have a too large mass. 
Moreover, the glueball should mainly decay to light hadrons, but not D*D*. Nevertheless, to 
measure the production angular distribution parameter a, where the differential cross section 
is proportional to (1 + a cos'^9), will be useful to clarify the glueball interpretation (with a 
negative value of a for the 0+"*' glueball). 

Could the X(4160) be an exotic charmonium-hybrid? say, the 1 hybrid, a possible partner 
of the 1 charmonium-hybrid ccg candidate, the Y(4260)[2Sl ED]- However, the problem is, 
a hybrid does not seem to have a favorable production mechanism associated with the J/ip 
in e^e~ annihilation. Compared with the production of double charmonium states such as 
e"*" + e" J/ip + ric{lS), the production of a ccg hybrid associated with J/ip requires an 
additional gluon production, and could therefore be relatively suppressed. But experimentally, 
the production rate of X(4160) J/-?/' is comparable to that of ric{lS)J/ip- 

Whether the X(4160) can be a charmonium molecule? The well know charmonium-like 
state X(3872) has been suggested being a D^D*^ + c.c. molecule either as a real bound state 
(see, e.g. [31] and references therein) or as a virtual state (see, e.g. [32|). The most significant 
motivation for the molecule assignment is that the mass of X(3872) is very close to the D^D*^ 
threshold. However, in the case of X(4160), its mass is above the D*D* threshold by about 
140 MeV. This makes the X(4160) unlikely to be a D*D* + c.c. molecule, since for molecules 
the binding energies due to meson exchanges are much less than 100 MeV. 

In summary, we have discussed various interpretations of the X(4160), observed by Belle in 
the process of double charm production e~^e~ — J/V'-|-X(4160) followed by X(4160) —>■ D*D*. 
The available information for this state from the data is its mass, width, a major decay mode, 
and its large production rate (comparable to ric{lS)) associated with J/ip in e'^e' annihilation. 

Using the leading order NRQCD calculation of the relative cross sections of double charmo- 
nium production in e^e~ annihilation as a guide, we find that though the D-wave spin-singlet 
charmonium state ^D2{2D) (with J^'~^ = 2~~^) could be a candidate of X(4160), the calculated 
production rate is too small, only about 5% of that for e+e~ J/ip + '7c(l'S'), and therefore 
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this assignment is unlikely. 

The possibility of X(4160) being the known '?/'(4160) produced via two-photon fragmentation 
in e^e~ annihilation is also discussed, but the calculated rate is much smaller (even with the 
3S-2D mixing effect included) than that for e^e~ — >■ J /ip + J /ip which, however, only has a 
small experimental upper limit. So the ?/'(4160) interpretation for X(4160) should be completely 
ruled out. 

The X(4160) is unlikely to be the excited 2"*"+ or l^^ charmonium state, e.g., Xc2(2P, 3P) or 
Xci(2P, 3P), because of the experimental non-observation and the calculated smallness of the 
cross sections for the J /'ip + Xc2(l-P) and J /i^ + Xci(l-P) production. 

The candidates of glueballs, ccg hybrids, and charmonium-molecules for the X(4160) might 
also be considered, but these interpretations are not very likely. 

In contrast to above interpretations, the ?7c(4S') assignment for X(4160) is an interesting 
possibility. The production rate of e+e~ J /ip + Vc{4:S) relative to e~^e~ ^ J/ip + Vci^S) 
in NRQCD is not very small, and could be compatible with the Belle data (note that for the 
observed e"^e^ ^ J/ip + ric{lS) cross section only a lower bound of 25.6 ± 2.8 ± 3.4 fb is given 
by Belle). And the non-observation of the DD mode of X(4160) can also be understood for this 
assignment. But one has to understand why ric{4:S) has such a low mass. And if one accepts 
X(4160) being the r]c{4:S), then the '?/'(4415) can hardly be the ip{4:S) as conventionally classified 
in charmonium spectrum. 

The Xco(3-P) is an even more interesting candidate for the X(4160). In particular, if the 
observed broad peak around 3.8-3.9 GeV in the recoil mass of DD against J/ip in + e~ —>■ 
J /ifj + D + D due to the Xco(2-P) state, then the Xco(3-P) assignment for X(4160) will be 
favored. However, as in the case of ric{^S) discussed above, one has to understand the problem 
of low mass values of the 3P states in the Xco(3-P) assignment, compared with conventional po- 
tential model calculations. We also emphasize that measurements on the angular distributions 
of cross sections are useful to distinguish between the Xco(3-P) and ?7c(4S') assignments for the 
X(4160). 

In order to clarify the nature of X(4160), it will be helpful experimentally to measure 
the differential cross sections (the production angular distributions), which can be different in 
different assignments, and to measure the strong decay branching ratios into various charmed 
meson pairs, and to measure the quantum numbers of X(4160). Theoretically, it is certainly 
important to have a reliable calculation for the strong decay rates, which is not very easy 
considering the complexity due to the coupled channel effects, and to understand why X(4160) 
has a dominant decay mode into D*D*. As for the production, as far as NRQCD is concerned, 
the NLO radiative corrections are only available for the e^e~ J/ip + r]c process[T5], which is 
confirmed by a recent independent calculation [33 j. It will certainly be very useful if the NLO 
calculation for the P-wave and D-wave charmonium states involved in double charmonium 
production can be performed. 

At present, we only have a limited understanding for the puzzling state X(4160). Since its 
finding was reported more than five months ago, there have been no theoretical papers on its 
interpretations. So, it is our hope that the discussion presented in this paper will stimulate 
more interesting discussions on this new hadronic state. 
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Table 1: Angular distributions and cross sections for double charmonium production in e+e~ 
annihilation at ^/s — 10.6 GeV with both QCD and QED contributions in the leading order 
NRQCD calculations (numbers without QED contribution are given with square brackets, see 
also text for the input parameters). 



Final state 


Differential cross section (fb) Cross section (fb) 


J/^ + r7e(15) 


2.47 [2.06]{1 + cos^e) 6.6 [5.5] 


J/i) +1 D2{ID) 


0.077 [0.069] (l + cos^^) 0.21 [0.19] 




0.123 [0.111] (l + cos2^) 0.34 [0.31] 


J/ij + r,,{AS) 


1.14 [0.95] (1 + cos2^) 3.0 [2.5] 




4.7 [4.6](l + 0.252cos2^) 10.2 [9.9 
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